Depletion forces between two spheres in a rod solution. 
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Abstract 

We study the depletion interaction between spherical particles of radius R 
immersed in a dilute solution of rigid rods of length L. The computed inter- 
action potential is, within numerical accuracy, exact for any value of L/R. In 
particular we find that for L ^ R, the depth of the depletion well is smaller 
than the prediction of the Derjaguin approximation. Our results bring new 
light into the discussion on the lack of phase separation in colloidal mixtures 
of spheres and rods. 



PACS: 82.70.Dd; 64.75.+g 
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Mixtures of colloids are abundant in industry as paints, glues, lubricants and other 
materials [|I|. They are also present in the preparation of foods and drugs, and in the 
biological realm: many living organisms or components of living organisms are colloidal 
suspensions of a variety of sizes and shapes. A pervasive issue in the colloidal domain is 
the stability of colloidal suspensions. Stability is often necessary for practical purposes as 
in paint formulation for instance, but certain applications like water treatment or mineral 
recovery might instead require aggregation or fiocculation to occur. The determination of 
the stability criteria or the study of the fiocculation kinetics can be achieved with reasonably 
good accuracy once the interparticle interaction potentials are known [Q. In a system of a 
pure colloidal species the two main interactions are the van-der-Waals attraction and the 
hard-core repulsion. Such a system is intrinsically unstable, the van-der-Waals attractive 
component of the potential always leading to fiocculation. In practice the stabilization of the 
suspension is enforced by using the screened electrostatic repulsion in aqueous solutions or 
steric polymer layers in organic solvents. To a good approximation the stabilized suspension 
can then be regarded as hard-core particles with no attractive potential component. One 
way of treating the stability issue in a colloidal mixture of two components is by considering 
the effective potential that the second species induces between two particles of the first 
species. 

A well studied case is the attractive potential that a solution of hard-core spheres 
of radius rg induce between two spheres of radius R, (see fig. ||): 

U,{H) = -kBT-<P-(2- + + -for H< 2r, (1) 

^ ' 8 ^ro V rj \ 3R 6RJ' " ^ ^ 

where H is the distance between the two large spheres and the volume fraction of the 
small spheres. Of course, Ug is zero if if > 2ro; by definition of the depletion potential the 
depletion at if = 2ro is the zero-point. We will only give the potential in the regions where 
it is non-zero from now on. The depletion attraction has been known since the pioneering 
work 1^ of Asakura and Oosawa who calculated the interaction energy between two flat 
plates immersed in a hard-sphere solution: 
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The depletion potential has a purely entropic origin: the exclusion of the small particles from 
the gap creates a pressure deficit and hence an effective attraction between the plates. When 
the exact analytical form of the interaction is only known for flat plates, one can still compute 
the interaction potential between spheres much larger than the range of the potential by 



Note that the Derjaguin approximation underestimates the attraction in this case. Indirect 
experimental evidence for the depletion effect were since long ago provided by the flocculation 
observed in colloids and emulsions, but more direct observations of the attraction have 
been only recently performed by force measurements [0 with surface force apparatus or 
microscopy-techniques . 

A much less-studied case is the depletion induced by hard-core rods, a surprising fact 
considering that rod-shaped particles in the colloidal range are present in a large variety 
of mineral and organic systems 0; they also exist in the biological realm where examples 
range from TMV-like virus to fibrils of amyloid /5-protein, the molecular agent at the origin 
of the Alzheimer disease. The depletion interaction between two flat surfaces immersed in 
a dilute rod solution was studied in |Q: 



with D the diameter of the rod, and 0* the Onsager volume fraction, below which the 
solution is isotropic. For spheres in a rod solution the potential can be written in general 
as: 



employing the so-called Derjaguin approximation 0. For instance, the asymptotic behaviour 
(as -R/(2ro) ^ 1) of eq. (|1]) can be obtained from eq. (^: 



UDerXH) = -7TR dH'u f{H') = -knT- (j) - {2 - -) (3) 
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Us{H) = kBT4.2 
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but no explicit form of Ki is available except in the Derjaguin approximation: Ki'^^'{x) = 
— (tt/G) (1 — x)^, which is expected to apply only for L/R <^ 1. Since this potential has 
quite a large contact value compared to ksT, one would expect on this basis to observe 
flocculation also in rod/sphere mixtures. Surprisingly, such a phase separation has not yet 
been clearly identified experimentally. A possible explanation for this lack of experimental 
evidence has recently been proposed [l^, based on repulsive contributions from the surface 



modification of the rod-rod excluded volume interactions. However such contributions are 



second order in the rod-concentration, and thus small below the Onsager concentration |]lT 
pI = A.2/{L'^D). The authors of reference |jlO| also note that most systems of interest involve 



rods that are of comparable size to the spherical particles; in particular the experiments |T2 
cited in [0] have 1.16 < L/R < 4.34, a range where the use of the Derjaguin approximation 
to compute the depletion potential in questionable. 

In this letter we present exact results for the depletion potential to first order in rod 
density. Our results are numerical in nature, but exact, i.e. the inaccuracy of our calculation 
is only limited by the numerical nature of the integration performed. It is for our case 
reducible to levels at which the error bars on the figures given below are not visible. The 
method of our calculation is based on the work on surface tension of objects immersed in 



rigid rod solutions that we published elsewhere ||T3|,|Tj]. We now briefly discuss this method 
and present our results. 

We consider two spheres of radius i? in a solution of rods of length L, and thickness D, 
with L/D ^ 1. We work to first order in the density of the rods, hence we really consider 
one rod and two spheres. A priori this is only valid when the rod-rod interactions can be 
ignored: for rod concentrations below the threshold pi- 

The interaction potential is given by the differences in the grand potentials of the 
spheres at infinite separation and inside the range of the potential, L. These grand potentials 
are integrals over the phase space available to a rigid rod: If the rod's center of mass is at 
a perpendicular distance z greater than L/2 to the surface of a sphere, then the rod is free 
to rotate, otherwise its rotational degrees of freedom are limited. The "depletion region" of 



the rod is shown in fig. |T]. This region can be further broken into two: In the first one (I) 
the rod interacts only with one sphere, whereas in the second one (II) interactions with both 
spheres are possible. Our previous work |jl3[ provides exact expressions for the contribution 
of region I to the free energy for all values of (L/R). Region II needs to be further broken 
down to three sub-regions where the functional form of the phase space allowed to the rod 
is different. The integration over the angular degrees of freedom can still be performed 
analytically, but integration over the possible center of mass positions is complicated by the 
lack of analytical expressions for the boundaries between the sub-regions where different 
functional forms apply. We therefore use the following procedure: once the position of the 
center of mass is given, we compute two angles, 9a and 9b, (see fig. ||), from its coordinates, 
and then numerically decide what functional form to use. Since we have exact expressions 
for the allowed phase spaces, this procedure reduces the problem to a two-dimensional 
numerical integration from a four-dimensional one, which can be performed rapidly on a 
personal computer. 

One can, of course, do the calculation using Monte-Carlo integration, and we have 
checked our results using this technique. When region II is not too small as compared 
to the total depletion region, Monte-Carlo yields reasonable accuracy, moderately fast. The 
results for these small H values are in perfect agreement, (see fig. with the exact method, 
for values of {L/R) anywhere from 0.1 to 15.0. When H gets larger, though, Monte-Carlo 
results fiuctuate considerably; longer machine time is needed to reduce these fiuctuations. 

Representative results for the function Ki{x) defined in eq. (|^) are shown in fig. ^ The 
Derjaguin function, K^^^'{x) is also shown for comparison. As can be seen from this plot, 
Derjaguin is indeed an excellent approximation to the depletion potential when {L/R) <^ 1. 
What is noteworthy is the large and systematic deviation from this approximation when 
{L/R) grows. This deviation is about 50% when {L/R) = 1, and becomes more pronounced 
as {L/R) increases. This can be understood in the light of general results we presented |T^ 
recently. These results relate the second derivative of the free energy of a rod system to 
the measure of configurations where the rod can wedge between two points on the surface. 



Clearly, in our case for a fixed H, this measure goes to zero as the rod gets longer. This 
implies that the free energy becomes a linear function of {L/R), and the function Ki, being 
related to the free energy by a factor of {L/RY, vanishes as R/L. Figure ^ shows the 
behaviour of the contact value, Ki{0), as a function of the rod/sphere ratio. We fitted the 
data of fifteen points to a function with three parameters. The resulting fit is: 

r.(r.._ ^ 1 + 0.8762 (L/i?) 

^ 6 1 + 1.33198 (L/i?) + 0.98225 (L/i?)2 ^ ^ 

as can be seen in fig. |^ the fit is very good indeed for the range 0.05 < L/R < 13. 



The lack of phase separation in the experiments [|12| is less surprising in the light of our 
results. Previous expectations based on Derjaguin approximation estimate the value of the 
attraction minimum at 10 — 20^5^- Our results show that the real value is much lower. This 
reduction amounts to a factor of three to five for these experiments and reduces the depth 
of the depletion potential to a few ksT. Note also that other factors like chain flexibility 
are likely to reduce the strength of the attraction further. 

We thank Phil Pincus for useful discussions. KY acknowledges support from the National 
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FIGURES 




FIG. 1. The depletion region: a) for H > L/2, b) for H < L/2. 
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FIG. 2. The basic construction of the numerical integration method. 
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FIG. 3. Monte-Carlo compared to the numerical integral for L/R = 0.1. 



10 




FIG. 4. The function Ki{x) for several values of L/R. 
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FIG. 5. Contact values.The function plotted here is the one given in the text. Notice the drastic 



change in the functional form of the data when L crosses 2R. 
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